Spin-orbit scattering in quantum diffusion of massive Dirac fermions 
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Effect of spin-orbit scattering on quantum diffusive transport of two-dimensional massive Dirac 
fermions is studied by the diagrammatic technique. The quantum diffusion of massive Dirac fermions 
can be viewed as a singlet Cooperon in the massless limit and a triplet Cooperon in the large-mass 
limit. The spin-orbit scattering behaves like random magnetic fields only to the triplet Cooperon, 
and suppresses the weak localization of Dirac fermions in the large-mass regime. This behavior 
suggests an experiment to detect the weak localization of bulk subbands in topological insulator 
thin films, in which a narrowing of the cusp of the negative magnetoconductivity is expected after 
doping heavy-element impurities. Finally, a detailed comparison between the conventional two- 
dimensional electrons and Dirac fermions is presented for impurities of orthogonal, symplectic, and 
unitary symmetries. 
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Recently, a featured observation in the quantum dif- 
fusive transport of topological insulators (TIs) [l|43| is 
the weak antilocalization (WAL) [4J-ll4|, an enhancement 
of conductivity due to destructive interference between 
time-reversed scattering loops. It exhibits as a negative 
cusp in low-field magnetoconductivity or logarithmically 
decreasing conductivity with increasing temperature in 
the metallic regime. 

The destructive interference is due to the ir Berry phase 
of massless Dirac fermions 



only the triplets due to spin relaxation, resulting in the 
crossover to WAL [26L l3lM33| . Therefore, the spin-orbit 
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If Dirac fermions ac- 
quire a mass, the tt Berry phase will be changed, lead- 
ing to the crossover from WAL to its opposite, the weak 
localization (WL) [17H2Q]. Because both their surface 
and bulk states are Dirac fermions [2l|, topological in- 
sulators provide a platform where rich quantum trans- 
port phenomena of Dirac fermions can be explored. The 
WL-WAL crossover has been observed [11] for the mag- 
netically doped surface states [22—24]. Suppressed WAL 
was also observed in thin films [9|, LL3| , probably due to 
the finite-size effect of the surface states [18|, [25| . More 
interestingly, the bulk states of topological insulators al- 
ways have the mass, which is nothing but the energy gap. 
Unlike conventional electrons, where spin-orbit coupling 
always leads to WAL [26I429J , the bulk states of topolog- 
ical insulator could give rise to WL due to their massive 
Dirac fermion nature [30]. 

In the language of "Cooperon" (the vertex of maxi- 
mally crossed diagrams that gives the quantum interfer- 
ence correction to conductivity), the conventional two- 
dimensional electrons have two spin- 1/2 bands crossing 
the Fermi surface, giving four "Cooperon" channels: one 
"singlet" and three "triplets", depending on the total 
spin angular momentum of incoming and outgoing elec- 
trons. Each channel contributes a conductivity correc- 
tion of the same magnitude but the singlet (triplet) gives 
WAL (WL). For only elastic scattering, the triplets out- 
number the singlet by 2, leading to WL in total. On 
the other hand, strong spin-orbit scattering can suppress 



scattering may also suggest an approach to probe the 
WL of the bulk states of TIs. However, it is still unclear 
what will happen to massive Dirac fermions if spin-orbit 
scattering is introduced. 

TABLE I: "Cooperon" channels for the conventional elec- 
trons and Dirac fermions. Triplet (singlet) channel gives WL 
(WAL). Spin-orbit scattering only quenches the triplet chan- 
nels, leading to the crossover from WL to WAL for conven- 
tional electrons and the suppression of WL in the large-mass 
limit of Dirac fermions. 7 is the Berry phase. 
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In this work, we study the effects of spin-orbit scatter- 
ing on the quantum diffusive transport of massive Dirac 
fermions. For a single cone of Dirac fermions, the trans- 
port can be understood as a "singlet Cooperon" in the 
massless limit and a "triplet Cooperon" in the large-mass 
limit. The spin-orbit scattering has no effect on the "sin- 
glet", but behaves like random magnetic impurities to 
the "triplet". As a result, the spin-orbit scattering can 
drastically suppress the WL in the large- mass regime, 
while has little impact on the WAL in the small-mass 
regime. We expect a narrowing of the cusp of the nega- 
tive magnetoconductivity of WAL, by doping impurities 
of heavy non-magnetic elements (e.g., Au) onto thin films 
of topological insulators. Finally, a systematic compari- 
son is presented for the conventional and Dirac fermions 
in the presence of scattering by ordinary, spin-orbit, and 
magnetic impurities. 

Both the two-dimensional bulk and surface states of 



topological insulator can be described by the Hamilto- 
nian of massive Dirac fermions, 



H = d ■ a, 



(i) 



where d = (hvk y ,—hvk x ,d z ), a is the vector of Pauli 
matrices, ft is Planck's constant over 27r, and v is the 
effective velocity. d z = A/2 — Bk 2 gives the mass for the 
Dirac cones in the effective model of the two-dimensional 
bulk subbands [30| or the surface bands [34J| , and d z = 
A/2 for magnetically doped surface states [22I424J. The 
spin-orbit scattering is described by 



Uso(r) 



ft 
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3 ■ Vu(r - RD x p, (2) 



the magnetic scattering time l/r m = 2/r x + l/r z and the 
spin-orbit scattering time l/r so = 2/r SOiX + l/r SOjZ . 

The new findings can be presented by reviewing the 
previous results. The massive Dirac model in Eq. (pp) 
carries a Berry phase given by 7 = 7r(l — cos#), which 
can be tuned from ir to as cos# changes from to 1, 
corresponding to the massless and large- mass limits, re- 
spectively. If there were only the elastic scattering, a sin- 
gle massless (7 = n) Dirac cone will exhibit WAL in the 
quantum diffusive transport. By introducing the mass 
term (7 7^ 7r), which changes the interference scenario 
from destructive to constructive due to the change of tt 
Berry phase to 0,a crossover from WAL to WL will be 
expected [H 0, [HI . 



where u(r — Rf°) represents the random potential by an 
impurity located at Rf°. The strength of the spin-orbit 
scattering will be characterized by a length £ so . Shorter 
4o means stronger spin-orbit scattering. £ so = \/Dr so , 
where D is the diffusion constant, r so is the spin-orbit 
scattering time with l/r so = 2/t so , x + l/r so ,z, 

1 7rN f 



ft 



T~so,z 

1 TTNfi 



(l + cos 2 «)n so < z (kxk')2, 
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where Np is the density of states at the Fermi energy, n so 
is the concentration of impurities that induce spin-orbit 
scattering. u so ^ x / z are defined from u(r — Rf°) by assum- 
ing delta potentials, cos# = dpj ' \/d 2 F + (ft/u/ci?) 2 , <^f is 
the value of d z at the Fermi surface, Uf is the Fermi wave 
vector. In above definitions, we have replaced (k x k ) ■ 
by its average over momentum directions (k x k') 2 . This 
replacement is valid as long as £ so ^> £ e [35] , and £ e is the 
length that characterizes the elastic scattering, which is 
always the strongest scattering mechanism in the quan- 
tum diffusive transport. Although intrinsic spin-orbit 
coupling carried by band dispersions of Dirac fermions 
and extrinsic spin-orbit scattering share the same ori- 
gin at the atomic level, their effects on the spin relax- 
ation 36|, [37| are quite different .similar to the difference 
between the Elliott- Yafet [38|, [39| and Dyakonov-Perel 
mechanisms I40|. 



With the help of the diagrammatic technique [26|, 1271 



41|-[44|, we calculated the quantum interference correc- 
tion to conductivity and magnetoconductivity for the 
massive Dirac model in the presence of elastic, spin- 
orbit, and magnetic scatterings. The formulas in the 
presence of spin-orbit scattering remain the same two- 
term structure as the one with only elastic and mag- 
netic scatterings [19|. The difference is the replacements 
1/7-3 -► I/t^-I/tso,^ and l/r z -> l/r z -l/r SOyZ , where r x 
and r z are the in-plane and out-of-plane magnetic scat- 
tering times, respectively. Besides, the total relaxation 
time now is given by 1/r = l/r e + l/r m + l/r so , where 
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FIG. 1: (a) Zero- field quantum interference correction to con- 
ductivity cr F (0) (in units of e 2 /ft) as a function of 7 and the 
spin-orbit scattering length l so . The dashed curve separates 
the positive WAL and negative WL regimes, (b) Magneto- 
conductivity Aa(B) for different 7 at £ so = 300nm. Other 
parameters: elastic scattering length ^ e =10nm, phase coher- 
ence length £(/) = 300 nm, and magnetic scattering length 

£m -> OO. 

We find that the spin-orbit scattering will suppress the 
weak localization in the large mass regime. Fig. HJb) 
presents the magnetoconductivity. WAL (WL) can be 
recognized as a sharp negative (positive) low-field mag- 
netoconductivity cusp. By introducing a small spin-orbit 
scattering (^ so = 300 nm), the negative curves of WAL 
with j/tt =1, 0.9, 0.75 are almost unaffected while the 
positive curves with 7 =0.5, 0.1, 0.001 are suppressed 
drastically. Even, the magnetoconductivity of WL to- 
tally vanishes in the 7—^0 limit. In Fig. [TJa), we present 
the quantum interference conductivity correction cr F (0) 
as a function of 7 and £ so . WAL (WL) demonstrates as 
positive (negative) cr F (0) at low temperatures and zero 
field. In the absence of spin-orbit scattering (^ so — >• 00), 
<j F (0) reaches positive maximum (WAL) when 7 = tt 
and negative minimum (WL) when 7 = 0. The bound- 
ary (dashed curve) that separates the positive WAL and 
negative WL regimes approaches to a saturate value of 
7 ~ 0.3 [T3] as £ so — »■ 00. As we increase the spin-orbit 



scattering (decreasing 4 ), the WL regime shrinks. In 
the limit of strong spin-orbit scattering (£ so — >> 0), the 
regime of WL vanishes, and the WAL-WL crossover on 
the £ so — »> oo side is replaced by WAL . 



The physical picture of the suppression of WL can be 

. The quantum 
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FIG. 2: (a) The weak antilocalization phase coherence length 
(^al) and the weak localization phase coherence length (^l) 
as functions of 7 for weak (£ so = lOOOOnm) and strong (£ so = 
300nm) spin-orbit scattering, (b) The weak localization phase 
coherence length ^l as a function of the spin-orbit scattering 
length £ so for different 7. 

The suppression of WL by spin-orbit scattering can 
be further studied by examining the effective phase co- 
herence lengthes. In the quantum diffusion transport, an 
electron can be scattered by static centers for many times 
but still maintains its phase. The phase, is protected 
by a long phase coherence length £$ (set as 300 nm in 
this work according to experimental fittings). In our cal- 
culated quantum interference correction to conductivity, 
the WL-WAL crossover can always be described by two 
competing terms, one for WL and the other for WAL. 
Each term looks like a Hikami-Larkin-Nagaoka formula 
[27j , and is characterized by an effective phase coherence 
length, denoted as £\, and £al, respectively. For very long 
£ so (weak spin-orbit scattering), £al diverges as 7 —> n, 
and £l diverges as 7 — >> [see dashed curves in Fig. |2fa)]. 
The divergence of effective phase coherence lengthes pro- 
tects the WAL as 7 — )> ir and WL as 7 — » 0. For short £ so 
(strong spin-orbit scattering), the WAL phase coherence 
length ^al remains almost unchanged in the small mass 
regime. However, the divergence of £\, in the large mass 
regime is suppressed [the solid curves in Fig. |2fa)]. Fig. 
EJb) shows £\, as a function of £ so in the large mass limit. 
The slope of the curves shows that the spin-orbit scat- 
tering has stronger influence on larger mass cases. In the 
7^0 limit, £l —> as £ so — >• 0, leading to the vanishing 
magnetoconductivity in Fig. QJb). 



understood by the "Cooperons" 
diffusion transport of a single cone of Dirac fermions can 
be understood as a "singlet Cooperon" in the massless 
limit and a "triplet Cooperon" in the large-mass limit, 
they give the two-term formulas for the quantum correc- 
tion to conductivity [19|. The spin-orbit scattering can 
be seen as spin-dependent magnetic fields. It has no effect 
on the singlet, because singlet state has zero total spin 
angular momentum of incoming and outgoing scattered 
electrons. However, spin-orbit scattering behaves like 
random magnetic fields to the triplet Cooperon, which 
carries a net total spin angular momentum of 1. In 
this way, the phase coherence length of large-mass Dirac 
fermions is shortened by spin-orbit scattering, much like 
by magnetic impurities. 
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FIG. 3: (a) A topological insulator thin film doped with impu- 
rities (e.g., Au) that bring strong spin-orbit scatterings, (b) 
The expected narrowing of the negative magnetoconductivity 
cusp after Au impurities are doped. In the simulation, one 
surface band (7 = it) and two bulk subbands (7 = 0.157r) are 
assumed. 



The suppression of WL by spin-orbit scattering pro- 
vides an experimental signature to detect the WL of bulk 
states in topological insulators. It is known that both the 
bulk and surface states contribute to the transport in as- 
grown topological insulators [45|, |46| . The experimentally 
observed negative magnetoconductivity was suggested to 
be the summation of the WAL from the surface and WL 
from the two-dimensional bulk subbands [30|. By intro- 
ducing impurities that can bring spin-orbit scattering, 
WL from the large-mass bulk states will be suppressed, 
but WAL from massless surface states will not be af- 
fected. Overall speaking, one will expect a narrowing of 
negative magnetoconductance cusp as shown in Fig. 03 
The recent experiment on the five quintuple-layer Bi2Te3 
thin film with Au impurities may already imply the effect 

0. 



TABLE II: Two-dimensional quantum diffusive transport of conventional and Dirac fermions for impurities of orthogonal 
(elastic), unitary (magnetic), and symplectic (spin-orbit) symmetries 47]. r m and r so are magnetic and spin-orbit scattering 
times, respectively. Elastic scattering (l/r e 7^ 0) is present in all the cases. 7 is the Berry phase. 



Conventional electron 

Massless Dirac fermion 7 = ix 

Massive Dirac fermion 7 6 (0, tt) 

Dirac fermion in large- mass limit 7 = 



Orthogonal 


Unitary 


Symplectic 


(l/r m = 1/t S o = 0) 


(l/r m / 0, 1/tso = 0) 


(l/Tm = 0, 1/Tso / 0) 


WL 


both suppressed 


WL-WAL crossover 


WAL 


suppressed WAL 


WAL f 


WAL-WL crossover 


both suppressed 


WAL-WL crossover^ 


WL 


suppressed WL 


suppressed WL* 



To summarize, in Table. Owe compare the quantum 
diffusive transport of 2D massive Dirac fermions with the 
conventional electrons in the presence of impurities of or- 
thogonal (elastic scattering), symplectic (spin-orbit scat- 
tering), and unitary (magnetic scattering) symmetries. 
The new results in this work are marked by f signs. 

This work is supported by the Research Grant Council 
of Hong Kong under Grant No. HKU 7051/11P. 

* Corresponding author: luhz@hku.hk 

t Corresponding author: sshen@hkii.hk 



[1] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

(2010). 
[2] X. L. Qi and S. C. Zhang, Rev. Mod. Phys. 83, 1057 

(2011). 
[3] J. E. Moore, Nature 464, 194 (2010). 
[4] J. G. Checkelsky, Y. S. Hor, M. H. Liu, D. X. Qu, R. 

J. Cava, and N. P. Ong, Phys. Rev. Lett. 103, 246601 

(2009). 
[5] H. L. Peng, K. J. Lai, D. S. Kong, S. Meister, Y. L. 

Chen, X. L. Qi, S. C. Zhang, Z. X. Shen, and Y. Cui, 

Nat. Mater. 9, 225 (2010). 
[6] J. Chen, H. J. Qin, F. Yang, J. Liu, T. Guan, F. M. Qu, 

G. H. Zhang, J. R. Shi, X. C. Xie, C. L. Yang, K. H. 

Wu, Y. Q. Li, and L. Lu, Phys. Rev. Lett. 105, 176602 

(2010). 
[7] J. G. Checkelsky, Y. S. Hor, R. J. Cava, and N. P. Ong, 

Phys. Rev. Lett. 106, 196801 (2011). 
[8] H. T. He, G. Wang, T. Zhang, I. K. Sou, G. K. L. Wong, 

J. N. Wang, H. Z. Lu, S. Q. Shen, and F. C. Zhang, Phys. 

Rev. Lett. 106, 166805 (2011). 
[9] M. H. Liu, C. Z. Chang, Z. C. Zhang, Y. Zhang, W. Ruan, 

K. He, L. L. Wang, X. Chen, J. F. Jia, S. C. Zhang, Q. 

K. Xue, X. C. Ma, and Y. Y. Wang, Phys. Rev. B 83, 

165440 (2011). 
[10] J. Wang, A. M. DaSilva, C. Z. Chang, K. He, J. K. Jain, 

N. Samarth, X. C. Ma, Q. K. Xue, and M. H. W. Chan, 

Phys. Rev. B 83, 245438 (2011). 
[11] M. H. Liu, J. S. Zhang, C. Z. Chang, Z. C. Zhang, X. 

Feng, K. Li, K. He, L. L. Wang, X. Chen, X. Dai, Z. 

Fang, Q. K. Xue, X. C. Ma, and Y. Y. Wang, Phys. Rev. 

Lett. 108, 036805 (2012). 
[12] J. Chen, X. Y. He, K. H. Wu, Z. Q. Ji, L. Lu, J. R. Shi, 

J. H. Smet, and Y. Q. Li, Phys. Rev. B 83, 241304(R) 



(2011). 
[13] Y. S. Kim, M. Brahlek, N. Bansal, E. Edrey, G. A. 

Kapilevich, K. Iida, M. Tanimura, Y. Horibe, S. W. 

Cheong, and S. Oh, Phys. Rev. B 84, 073109 (2011). 
[14] H. Steinberg, J. B. Laloe, V. Fatemi, J. S. Moodera, and 

P. Jarillo-Herrero, Phys. Rev. B 84, 233101 (2011). 
[15] T. Ando, T. Nakanishi, and R. Saito, J. Phys. Soc. Jpn. 

67, 2857 (1998). 
[16] D. Hsieh, Y. Xia, D. Qian, L. Wray, J. H. Dil, F. Meier, 

J. Osterwalder, L. Patthey, J. G. Checkelsky, N. P. Ong, 

A. V. Fedorov, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, 

R. J. Cava, and M. Z. Hasan, Nature (London) 460, 1101 

(2009). 
[17] K. I. Imura, Y. Kuramoto, and K. Nomura, Phys. Rev. 

B 80, 085119 (2009). 
[18] P. Ghaemi, R. S. K. Mong, and J. E. Moore, Phys. Rev. 

Lett. 105, 166603 (2010). 
[19] H. Z. Lu, J. R. Shi, and S. Q. Shen, Phys. Rev. Lett. 

107, 076801 (2011). 
[20] V. K. Dugaev, P. Bruno, and J. Barnas, Phys. Rev. B 

64, 144423 (2001). 
[21] H. J. Zhang, C. X. Liu, X. L. Qi, X. Dai, Z. Fang, and 

S. C. Zhang, Nat. Phys. 5, 438 (2009). 
[22] Y. S. Hor, P. Roushan, H. Beidenkopf, J. Seo, D. Qu, J. 

G. Checkelsky, L. A. Wray, D. Hsieh, Y. Xia, S.-Y. Xu, 

D. Qian, M. Z. Hasan, N. P. Ong, A. Yazdani and R. J. 

Cava, Phys. Rev. B 81, 195203 (2010). 
[23] Y. L. Chen, J. H. Chu, J. G. Analytis, Z. K. Liu, K. 

Igarashi, H. H. Kuo, X. L. Qi, S. K. Mo, R. G. Moore, 

D. H. Lu, M. Hashimoto, T. Sasagawa, S. C. Zhang, I. 

R. Fisher, Z. Hussain, and Z. X. Shen, Science 329, 659 

(2010). 
[24] L. A. Wray, S. Y. Xu, Y. X. Xia, D. Hsieh, A. V. Fedorov, 

Y. S. Hor, R. J. Cava, A. Bansil, H. Lin, and M. Z. Hasan, 

Nat. Phys. 7, 32 (2011). 
[25] Y. Zhang, K. He, C. Z. Chang, C. L. Song, L. L. Wang, 

X. Chen, J. F. Jia, Z. Fang, X. Dai, W. Y. Shan, S. Q. 

Shen, Q. Niu, X. L. Qi, S. C. Zhang, X. C. Ma, and Q. 

K. Xue, Nat. Phys. 6, 584 (2010). 
[26] G. Bergmann, Phys. Rep. 107, 1 (1984). 
[27] S. Hikami, A. Larkin,and Y. Nagaoka, Prog. Theor. Phys. 

63, 707 (1980). 
[28] S. Maekawa and H. Fukuyama, J. Phys. Soc. Jpn. 50, 

2516 (1981). 
[29] G. Bergmann, Phys. Rev. Lett. 48, 1046 (1982). 
[30] H. Z. Lu and S. Q. Shen, Phys. Rev. B 84, 125138 (2011). 
[31] A. J. Millis and P. A. Lee, Phys. Rev. B 30, 6170 (1984). 
[32] K. B. Efetov, A. I. Larkin, and D. E. Khemlnitzkii, JETP 



52, 568 (1980). 
[33] B. L. Altshuler, A. G. Aronov, A. I. Larkin, and D. E. 

Khemlnitzkii, JETP 54, 411 (1981). 
[34] H. Z. Lu, W. Y. Shan, W. Yao, Q. Niu, and S. Q. Shen, 

Phys. Rev. B 81, 115407 (2010). 
[35] J. Rammer, Quantum Transport Theory (West view 

Press, Boulder, CO, 2004). 
[36] P. D. Dresselhaus, C. M. A. Papavassiliou, and R. G. 

Wheeler, Phys. Rev. Lett. 68, 106 (1992). 
[37] S. V. Iordanskii, Yu. B. Lyanda-Geller, and G. E. Pikus, 

JETP Lett. 60, 206 (1994). 
[38] R. J. Elliott, Phys. Rev. 96, 266 (1954). 
[39] Y. Yafet, Phys. Rev. 85, 478 (1952). 
[40] M. Dyakonov and V. Perel, JETP 33, 1053 (1971). 
[41] H. Suzuura and T. Ando, Phys. Rev. Lett. 89, 266603 



(2002). 
[42] E. McCann, K. Kechedzhi, V. I. Fal'ko, H. Suzuura, T. 

Ando, and B. L. Altshuler, Phys. Rev. Lett. 97, 146805 

(2006). 
[43] B. L. Altshuler, D. KhmePnitzkii, A. I. Larkin, and P. A. 

Lee, Phys. Rev. B 22, 5142 (1980). 
[44] N. H. Shon and T. Ando, J. Phys. Soc. Jpn. 67, 2421 

(1998). 
[45] Y. S. Hor, A. Richardella, P. Roushan, Y. Xia, J. G. 

Checkelsky, A. Yazdani, M. Z. Hasan, N. P. Ong, and R. 

J. Cava, Phys. Rev. B 79, 195208 (2009). 
[46] D. X. Qu, Y. S. Hor, J. Xiong, R. J. Cava, and N. P. 

Ong, Science 329, 821 (2010). 
[47] F. J. Dyson, J. Math. Phys. (N.Y.) 3, 140 (1962). 



